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We develop a composite Dirac fermion theory for the fractional quantum Hall effects (QHE) 
near charge neutrality in graphene. We show that the interactions between the composite Dirac 
fermions lead to dynamical mass generation through exciton condensation. The four-fold spin-valley 
degeneracy is fully lifted due to the mass generation and the exchange effects such that the odd- 
denominator fractional QHE observed in the vicinity of charge neutrality can be understood in terms 
of the integer QHE of the composite Dirac fermions. At the filling factor v = 1/2, we show that the 
massive composite Dirac fermion liquid is unstable against chiral p-wave pairing for weak Coulomb 
interactions and the ground state is a paired nonabelian state described by the Moore- Read Pfaffian 
in the long wavelength limit. 
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Since the discovery of graphene [T] , a rich set of inte- 
ger [2HS] and fractional P-P quantum Hall effect (QHE) 
have been explored experimentally in the new setting of 
a two-dimensional electron gas with relativistic disper- 
sion. Currently, the v = \/Z state has been convincingly 
observed [BHn]- Remarkably, initial evidence for the emer- 
gence oi & V — \/2 state has been observed [S], raising 
the hope for realizing nonabelian statistics in suspended 
graphene. 

For graphene, it is known that there is a spin-valley 
SU(4) symmetry. This restricts the filling factors v of the 
integer QHE of the Dirac fermions to ±2, ±6, ±10, • • • , 
provided that the SU(4) symmetry is unbroken. The 
observation of the v = ±.1 integer Hall plateaus [Sj in- 
dicates that the SU(4) symmetry is broken; dynamical 
mass generation, which lifts the spin-valley symmetry 
of the zeroth Dirac Landau level {(LL) [TU], as well 
as quantum Hall ferromagnetism |11H13) have been pro- 
posed for its explanation. Despite recent theoretical ef- 
forts on the fractional QHE (FQHE) in graphene [BHTB] . 
it remains uncertain whether the observed v = 1/3 state 
implies breaking of the SU(4) symmetry and is thus a 
single-component Laughlin state, or a multi-component 
Halperin state with spin- valley degeneracies '18 . More- 
over, although states at v = ±3/2, ±7/4 have been inves- 
tigated previously |19j. whether paired, nonabelian quan- 
tum Hall states can emerge at = ±1/2 of the (LL in 
graphene has not been studied. 

In this paper, we propose a mechanism of dynamical 
mass generation by exciton condensation and exchange 
driven polarization to describe both the abelian and the 
nonabelian FQHE in graphene. To this end, we extend 
the composite fermion Chern-Simons (CS) theory [201 
to the case of Dirac particles attached to an even number 
of flux quanta through the CS gauge field. A crucial 
step in our theory is to work with the proper particle 



density via a particle-hole transformation such that the 
vacuum state of the relativistic composite Dirac fermions 
(CDF) is defined by the charge neutral state with all 
negative energy states filled. The CDF theory introduced 
below are rather general, involving only relations between 
the CDF particle density and its filling fraction V which 
corresponds to a unique filling fraction of the original 
electrons once the ground state is determined. 

At the filling fractions i/ — with an even inte- 

ger, there are two different types of CS transformations: 
one preserves and one breaks the SU(4) symmetry. By a 
variational calculation, we show that it is energetically fa- 
vorable to develop an exciton condensate. The latter sup- 
ports single quasiparticle excitations with a CDF mass 
gap at low density above the condensate. If the mass 
gaps have the same sign for all SU(4) components, the 
SU(4) symmetry is preserved. This kind of CDF systems 
have a non-zero Chern number ±2, implying that the ex- 
citon condensate is an anomalous Hall liquid. In terms 
of the Bloch electrons in graphene, this corresponds to 
a total filling factor i' — ±(2 — ?). However, when the 
exciton mass gaps for different spin (or valley) compo- 
nents are of different signs, the SU(4) symmetry is re- 
duced to SU(2)(X)SU(2), i.e., the spin (or valley) degener- 
acy is lifted. In this case, the exchange term in the sta- 
tistical interaction further breaks the remaining SU(2) 
symmetry so that the ground state is fully spin-valley 
polarized and the exciton condensate becomes topolog- 
ically trivial. Such states have v — ±17. Thus, in our 
theory, the dynamical, SU(4) symmetry-breaking mass 
generation offers a route to the observed FQHE states in 
graphene near charge neutrality [BHS] • We show that the 
quasiparticles above the exciton condensate have a pair- 
ing instability in the chiral p-wave channel, leading to an 
even denominator paired quantum Hall state described 
by the Moore-Read Pffafian [551 US] in the long wave- 
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length limit. On the other hand, at odd-denominator 
filling fractions such as v = 1/3, the quasiparticles above 
the exciton condensate occupy fully filled Landau levels 
of the residual magnetic field [20 , giving rise to a single- 
component Laughlin state for the electrons in graphene. 

We start with graphene electrons in the con- 
tinuum limit under a perpendicular external mag- 
netic field. The effective Hamiltonian reads H — 
—ihvp J xtplai{di + ieAi/h^jtl^g where s is spin in- 
dex and Q!j = 7°7\ The 7 matrices are given by 
7O = / (g) n, 7' = — iiJi (8) T2 where ai and Ti are the 
Pauli matrices in the sublattice and valley subspaces, re- 
spectively. The CDF field 'i's is introduced by a uni- 
tary transformation V's(x) = e''^=^^^4's(x) with Is(x) = 
/ d^x'ICss'Ps' (x') arg(x— x') and Ps the particle den- 
sity. In terms of "^fg, H = ^g—ihvF J d^x'^laiidi + 
ie{Ai + af where the CS gauge field af is defined 

by V X a''(x) = ICgs' Ps'{^)h/e . The condition for 
CDF to experience a vanishing net magnetic field on aver- 
age requires K.^ (pi) + /C^^ (p^) = JC^i (pi) + K.^^ (p^) = 
eB/h. Among the possible solutions of the /C-matrix, 
two simplest and physically relevant ones at the filling 
fraction i' = {{pi) + {p^))h/eB = are given by 



\sgn(i?)0 -(j) J 



= sgn(?) 




(1) 

An obvious distinction between K} and /C^ is that the 
former breaks the spin rotation symmetry while the lat- 
ter preserves the full SU(4) symmetry. The solution as- 
sociated with K} requires (p^) = eB/h(f), {p^) ~ for 
u = 1/0 and (p^) = —eB/h(j), (pi) = for j/ = —1/0. 
Thus K.^ describes a state in which spin must be fully 
polarized. Equivalently, one can replace the spin with 
valley to break the valley symmetry. 

The CDF transformation introduces a statistical inter- 
action, which takes an instantaneous form, 

V,t=J2~'f'''F J d'a;dV (*l,,(x)A^,(x,x')*ii«(x) 



-*LW>(s(x,x')*L.(x: 

where R and L label the two valleys and 

l/{z-z') 
l/{z-z') 



(2) 
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i(x,x') ^^JCss'Sps'{x.') 



with Sps = ps ~ (ps) and the holomorphic coordinate 
z = x + iy. Since the CDF experiences zero net magnetic 
field, its field operator can be expanded in the helicity 
basis as 
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where = arctan /Cy/fca,. In this basis, the ki- 
netic energy part of the Hamiltonian becomes H = 
E.k fi«Ffc(AL^.k + i?Li3sk + Cl^Csv + Dl^Ds^). 

We next show that the statistical interaction drives the 
formation of an exciton condensate for the CDFs. The 
variational wavefunction of the exciton vacuum can be 
written as follows, 



s k 



|0> = ]J ]J (cos ipRsk - sin-^j^sk^L-^L 

' cos ipLs-k - sin ^LskC'l^Dl^) |vac) , (3) 



where fB.,L are the variational parameters and |vac) de- 
scribes the state where the valence CDF bands are filled 
and the conduction bands empty. The quasiparticle op- 
erators associated with |0) are given by, 
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sin(y5fl,jkSsk, 

-I- cos (pBskBsli, 



(4) 



for the i?- valley and similarly in terms of c^^. and dgk for 
the L-valley. The normal state with a nonzero particle 
density can be written as |iV) = Uk<k§ 4 Uk<k^ 4 |0) 
where kp' are the Fermi wave vectors. Hereafter, we 
measure the wave vector in unit of the spin-polarized 
CDF Fermi vector kp = 1/Ib4>^^^ and energy in hvpkp. 
We will show that the variational energy is indeed min- 
imized at nonzero (Pr.l- It is useful to note that the 
exciton condensation energy from the statistical interac- 
tion is linear in the exciton mass. This is essentially due 
to the fact that Vgt is a current-density interaction and 
can be expressed as. 



s,s' 

X / d\d\'4'3(x)7V*s'(x') 



2^j2 /,f, / \„,3„,5,T, f,J\(^ x') • Js's(x',x) 



(5) 

where Ji^s'<;(x', x) = 'l's'(x')7*^'s(x). The internal 
energy coming from the statistical interaction (Vst) 
is a linear function of the exciton order parameter 
(^'s(x)7^7^5's(x')). As a result, the mass has the same 
sign as that of the corresponding matrix elements in /C 
in order to gain energy. To gain further insight, one can 
isolate the interaction Vst in the exchange channel (e.g., 
for the _R-valley) 



-T 

s,k,p 



sm 2ipiisu^ cos 2(pRsk + t 
i^ss r flsk^skaipflsp 

(6) 



where k> (k<) is the bigger (smaller) of k and p. To 
lower energy, i.e. to have a negative F^t, JCgs and tpRsk 
must have the same sign and a positive tp implies a pos- 
itive mass. The diagonal elements of JC^ have opposite 
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signs, leading to opposite masses for the spin-down and 
spin-up CDF bands. In contrast, they have the same sign 
in /C^; and the mass gaps are all positive as illustrated in 
Fig.[T] 

Thus, the exciton condensation is essential to classify 
the different states described by the two types of K. ma- 
trices. The difference lies in the topology characterized 
by the Chern number of the CDF bands, which is de- 
termined by the sign of the mass. From the sign of the 
mass gap, one can read out the Chern number of differ- 
ent states described by the JC matrices: For K} , the oc- 
cupied bands has a zero total Chern number and so the 
total filling factor is given by v = v = ±l/(/) following 
the SU(4) symmetry breaking. For the SU(4) symmetric 
case described by K? , the identical mass gaps contribute 
to a non-vanishing total Chern number ±2, resulting in 
a quantum anomalous Hall effect (QAHE) for the CDF. 
Remarkably, this QAHE at v implies that the total filling 
factor V — ±(2 — 1/0) for the electrons, where the integer 
contribution to a^y comes from the QAHE of the exciton 
condensate. In the rest of the paper, however, our focus 
will be on the FQHE in the vicinity of charge neutrality 
observed in recent experiments [BHS]- 
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FIG. 1: (color online) The difference between K} and K? is 
reflected in the signs of the exciton mass gaps. K} leads to 
opposite masses for the two spin projections and therefore 
describes the u — state (top panel). K.^ results in the 
same sign of mass for all bands, implying filling fractions u = 
— 2 -|- 1/0 (bottom panel). 



We now calculate the dynamical mass for the sym- 
metry breaking state \N) described by /C^ at ly = 1/2 
by minimizing the variational energy i?jv = {N\H\N) 
where the Hamiltonian includes both the statistical in- 
teraction Vst and the Coulomb interaction V^(x,x') ~ 
g/jx — x'l {g — /ATrehvp)- We focus on the spin-down 
bands and drop the spin indices. The variational equa- 
tion for the i?-valley (puu can be expressed as a set of 
self-consistent equations for the quasiparticle dispersion 
eflk = \/«flk + "^Ifk' where rriRk is the mass gap, Q;_Rk 
is the renormalized dispersion, and sin2(pjjk = Wiik/e/fk 



(same for the L- valley), 

^ p ^«p ^p>feP^«p 

Here, — ^ — n^, is the occupation number of 
the i?- valley, and vi{k,p) is the coefficient of the angu- 
lar expansion of the Coulomb interaction Vc in the ^-th 
angular momentum channel. Since projects out filled 
states, only states above the Fermi level contribute to the 
dynamical mass. Note that a natural ultraviolet energy 
cutoff for ([?]) is the energy spacing between (LL and the 
first LL, which is the largest energy scale in the prob- 
lem. Restoring the unit, the mass can be expressed as 
MfOi — mi^if^hkp /vp and scales according to Mijk (x ^/B. 
It is worth emphasizing that, in contrast to the nonrel- 
ativistic composite fermion theory where an appropri- 
ate composite fermion mass remained elusive |21j . the 
mass of the relativistic CDF theory naturally arises from 
the interactions through exciton condensation. We found 
that, due to the exchange interaction, the valley polar- 
ized state {kp = \f2kp, kp = 0) has lower energy than 
the unpolarized state {kp = kp ~ kp). For example, for 
Coulomb strength g = 0.3, (j) — 2 and a momentum cut- 
off A = 2fci?, the energy density of the polarized state is 
approximately —0.88 while that of the unpolarized state 
is —0.85. The solution for the complete dispersion of the 
exciton mass m/jk in the spin-valley polarized state is 
shown in Fig. [3^. The resulting CDF band structure is 
shown schematically in Fig. [2j The unoccupied spin-up 
bands have a mass with the same magnitude but oppo- 
site sign as the unoccupied spin-down band dictated by 
the diagonal elements of K.^ discussed before. It should 
be stressed that the mass gap of the empty CDF bands 
generated by the statistical interaction is large enough 
so that the chemical potential lies inside the gap, mak- 
ing the emergence of the spin-valley polarized state fully 
self-consistent with /C^. 

Before turning to the pairing instability, we comment 
on odd-denominator filling fractions. In the case of 

— ±1/3, these /C matrices with = 2 does not lead 
to the complete cancellation of the external magnetic 
field. Rather, the CDFs fill the Landau levels of the 
residual magnetic field B* = B — 2\p\h/e at an integer 
filling. These CDF Landau levels will develop mass gaps 
through the statistical interaction and Coulomb interac- 
tion via magnetic catalyst at integer filling factors [lOj . 
Under exchange splitting, only one CDF Landau level 
will be completely filled and exhibits an integer QHE. 
This corresponds to a FQHE state at v — ±1/3 of the 
electrons [20] described by the single-component Laugh- 
lin state. 

We now show that, at even denominator filling frac- 
tions, the spin-valley polarized composite massive Dirac 
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FIG. 2: (color online) Schematic plot of the spin-down CDF 
bands before and after exciton condensation (EC) and ex- 
change induced polarization. 



fermion liquid has a pairing instability where the quasi- 
particles on top of the exciton condensate form spin- 
triplet pairs in the chiral p-wave channel. In terms of 
the i?-valley, the statistical pairing interaction is domi- 
nated by the 1=1 angular momentum channel and has 
the form, 



TT sin 2y?ji^k> cos 2v;j;sk< 

V ^ fc> 



(8) 



It is remarkable that this pairing interaction is present 
only if there is an exciton condensate, i.e., when lyJ^sk 7^ 
0. Since the latter requires Landau level mixing, this im- 
plies that Landau level mixing is crucial for the pairing 
to occur. The variational wavefunction for the paired 
state has the BCS form jfi) = HkC^k + v\^alM}_^ |0), 
with jwkp + |wkP — 1, where |0) is the exciton vacuum 
defined in ([s]). Note that the variational wavefunction 
contains both the exciton and pairing order parameters 
which must be determined self-consistently by minimiz- 
ing the ground state energy. The self-consistent equa- 
tions for the dynamical mass have the same form as in 
([t]) with = |ukP, but unlike in the normal state, ex- 
citon pairs exist even for k < kp due to pairing. The 
variation of the energy with respect to and leads 
to the familiar BdG equation (dropping the valley index) , 



-EkWk = CkMk + ^k^k, 

£^kWk = -CkWk + ^kUk, 



where Ck = ek - /9k 
/3k = ff5Ep«o(fc,p)(l 



(9) 



and 



,) comes from the Coulomb 



exchange. The gap function Z\k is determined by the gap 
equation 



V ^ 
p 



UpVp 



r 1 TOk> 'ak< 
k> ek> ek< 



1 



-I vi[k,p} - -——vo(k,p) 



CkEp 



2 EkCp 



Fig. 3 displays the numerical solution of i?k and Z\k at 
(j) = 2 and g = 0.3, where the ground state is indeed a 



chiral p + ip paired state. Because the Coulomb inter- 
action is pair-breaking, the pairing gap Age = 2min(£'k) 
reduces with increasing gjmd vanishes at a critical value 
gc ~ 0.53 and 1.28 for </> = 2 and 4 respectively. For 
g > gc, the massive Dirac fermions form a stable Fermi 
liquid state. 
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FIG. 3: (color online) (a) The dynamical masses m_Bk and 
m^k in the chiral p-wave paired state. The corresponding ex- 
citon masses in the normal state are very close and drawn in 
dashed lines, (b) Solution of the BdG equation for the pairing 
gap function zik, the condensate amplitude /k = u'^vi^, and 
the momentum distribution function rik; (c) The CDF disper- 
sion in the normal state and that of quasiparticle excitations 
in the paired state; (d) Pairing wave functions (unnormalized) 
in real space. All results are for (j> = 2 and g = 0.3 with a 
momentum cutoff A = 2k f- Due to valley splitting, the Fermi 
vector is kp = \/2kF- 

To gain further insights into the paired state, we study 
the two-particle pairing wave function by projecting the 
BCS state to real space: g(xi-X2) = (0|^'(xi)*(x2)|fi). 
Because of the spinor structure, the upper and lower com- 
ponents are obtained separately. 



51 W = y XI (cosy'flk + sin(/?_Rk)^ e ^^^'-e^'^-^Sk, 

k 

52 (x) = y X! i^'^^'^R'^ - sin(/3iik)^ e''' ''5k, (10) 

k 

where gk = '^'k/uk. In the long wavelength limit, 
g\^ sa e^^^/k and m/jk ~ 'n^o is essentially a constant 



(c.f. Fig. [3]). We have, to quadratic order, (cos(^i^k 
sin(yS_Rk)^ w 2 - v'p^k'^/2ml and {cos ipm^ - sin(p;?k)^ ~ 
v'p^k"^ /2mQ where v'p is the renormalizcd Fermi velocity. 
Thus, the large distance behaviors of the pairing wave 
function arc 51(7") cx 1/z + c-y/2/7r cos(Kr — 7r/4)/zy^ 
and g2{r) oc —c^y2jTT cos{Kr — 7r/4)/zy^, where c and 
K are numerical constants. Except for the subleading 
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oscillatory contributions that decay faster at large dis- 
tances, the pairing wave function resides predominantly 
on one component of the spinor and has the form of the 
Moore- Read Pfaffian state in agreement with the numer- 
ical results shown in Fig. [3] 

In summary, we presented a composite relativistic 
fcrmion theory for the FQHE in graphene. We showed 
that the ground state has spontaneous spin- valley polar- 
ization, leading to a Laughlin-like state at v — 1/3 and 
a chiral j3-wave pairing state of the Moore-Read Pfaffian 
at 1/ = 1/2. A crucial prediction of the present theory 
is the dynamical mass generation for the Dirac fermions 
through the formation of an exciton condensate in the ze- 
roth Landau level that originates from Landau level mix- 
ing and facilitates the spin-valley polarization and the 
pairing interaction for the paired states at even domi- 
nator filling fractions. This CDF mass scaling with \fB 
should be detectable by scanning tunneling microscopy 
in high magnetic field. In a recent study [8] on suspended 
graphene, a plateau like feature near = 1/2 is observed 
in some but not all samples. Future experiments are de- 
sirable to explore whether a true single-component quan- 
tum Hall state emerges at = 1/2 that realizes non- 
abelian statistics in graphene. 

This work is supported in part by DOE grant DE- 
FG02-99ER45747 and the NSF of China through the 973 
program of MOST. 
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